Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 2 February 2008 (MN KTpX style file v2.2) 



A softer look at MCG-6-30-15 with XMM-Newton 



A. K. Turner 1 *, A. C. Fabian 1 , S. Vaughan 1 and J. C. Lee 2 

1 Institute of Astronomy, Madingley Road, Cambridge CB3 0HA 

2 Chandra Fellow, Massachusetts Institute of Technology, Center for Space Research, 77 Massachusetts Ave. NE80, Cambridge, MA 
02139, USA 



2 February 2008 



o 
o 

(N 



(N 
> 
oo 

m 
O 
m 

o 
6 



ABSTRACT 

We present the analysis and first results from the Reflection Grating Spectrometer (RGS) 
during the 320 ks XMM-Newton observation of the Seyfert 1 galaxy MCG-6-30-15. The 
spectrum is marked by a sharp drop in flux at 0.7 keV which has been interpreted by 
Branduardi-Ravmont et al. using RGS spectra from an earli er and sho rter observation as the 
blue wing of a strong relativistic OVIII emission line and bv lLee et al.l using a Chandra spec- 
trum as due to a dusty warm absorber. We find that the drop is well explained by the Fel 
L2,3 absorption edges and obtain reasonable fits over the 0.32-1.7 keV band using a multi- 
zone, dusty warm absorber model constructed using the photoionization code CLOUDY. Some 
residuals remain which could be due to emission from a relativistic disc, but at a much weaker 
level than from any simple model relying on relativistic emission lines alone. A model based 
on such emission lines can be made to fit if sufficient (warm) absorption is added, although 
the line strengths exceed those expected. In order to distinguish further whether the spectral 
shape is dominated by absorption or emission, we examined the difference spectrum between 
the highest and lowest flux states of the source. The EPIC pn data indicate that this is a power- 
law in the 3-10 keV band which, if extrapolated to lower energies, reveals the absorption 
function acting on the intrinsic spectrum, provided that any emission lines do not scale ex- 
actly with the continuum. We find that this function matches our dusty warm absorber model 
well if the power-law steepens below 2 keV. The soft X-ray spectrum is therefore dominated 
by absorption structures, with the equivalent width of any individual emission lines in the 
residuals being below about 30 eV. 

Key words: galaxies: active - galaxies: Seyfert: general - galaxies: individual: MCG-6-30- 
15 - X-ray: galaxies 



1 INTRODUCTION 

MCG-6-30-15 is one of the most studied Seyfert 1 galaxies in 
the X-ray band due to its high X-ray flux and relative close- 
ness (z = 0.00775, L 2 -iokoV ~ 5 x 10 42 ergs _1 ; Ho = 
70kms~ 1 Mpc -1 ). At energies above 2 keV the X-ray spectrum 
consists of a power-la w with V ~ 1.9, skewed iron Kq line and 
reflection hump (e.g. Fabian et al. 20o3). At softer energies much 
controversy still exist concerning the observed spectral features. 
Before the use of grating spectra, sharp drops in flux observed in 
soft X-ray spectra were interpreted as absorption edges due to OVII 
and OVIII in ionized gas within the AGN environment ("the warm 
absorber" Totani et all 1 9961 lRevnoldil997tlGeorge et alll998t) . 

However the increase in spectral resolution afforded by the 
use of grating spectrometers onboard XMM-Newton and Chan- 
dra has challenged the view that reprocessing of soft X-rays by 
photoionized gas is responsible for the observed spectral features. 
iBranduardi-Ravmont et alj 1200 ll) proposed that sharp drops seen 
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in the Reflection Grating Spectrometer (RGS) spectra from XMM- 
Newton observations of MCG-6-30-15 and Mrk 766 at ~ 0.39, 
~ 0.55 and ~ 0.7 keV were due to the blue wings of relativistically 
blurred emission lines of CVI, NVII and OVIII Lya, the blurring 
mechanism being the same as that used to explain the skewness of 
the Fe Ka line. The reason for adopting this view was that the po- 
sition of the sharp drop at ~ 0.7 keV, if caused by an OVII absorp- 
tion edge inferred an unrealistically high infall velocity of matter 
of ~ 16, 000 kms -1 . iLee et all 1200 ll) countered this suggestion 
with Chandra High Energy Transmission Grating Spectrometer 
(HETGS) observations in which the 16,000 krns" 1 discrepancy in 
position of the OVII edge was explained by the blending together of 
resonance absorption lines red-ward of the edge and by L-shell ab- 
sorption from neutral iron which is probably in the form of dust (the 
"dusty warm absorber " model; D WA). In a re-analysis of the orig- 
inal RGS data lSakoetai]l2003l) fitted both the DWA and blurred 
emission line models to the original RGS data. These workers argue 
that the models of Lee et al. 1 2001) over-predicted the flux below 
the passband of the HETG instrument and that the relativistic line 
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interpretation provides a better fit to the data than the DWA inter- 
pretation. 

Given the evidence for a relativistically-blurred iron line 
and reflection continuum in MCG-6-30-15, some relativistically- 
blurre d soft X-ray lines are also expected in the spec- 
trum iRo ss & Fa bianll993tlNavakshin. Kazanas. & Kallmarl2000t 
iRozanskaetalJI 2002) The strength of these lines is however pre- 
dicted to be much smaller, and the blue wing broader, than those 
claimed from the RGS spectra iBallantvne. Ross. & Fabiar|2002l) . 
The determination of the underlying soft X-ray continuum is of im- 
portance for understanding the X-ray irradiation of accretion discs. 

This paper presents the RGS d ata for the l ong ~ 320 ks ob- 
servation of MCG-6-30-15 lFabianetalj r2002). In particular we 
test to see if a DWA or relativistically-blurred emission line model 
can explain the overall spectral shape and see what can be inferred 
about the underlying continuum. The rest of the paper is organised 
as follows. Section[2|describes the data reduction. In Section|5]ab- 
sorption across the 0.7 keV drop is examined. Section |4] presents 
a photoionization model of the warm absorber used in Section [5] 
where DWA and the relativistically-blurred emission line models 
are fit to the data. In Section |6| an attempt is made to use spectral 
variability to determine the underlying continuum. The results are 
discussed in Section|7| This paper does not represent a full descrip- 
tion of the data but a first step in modelling the RGS spectrum of 
the DWA in MCG-6-30-15 and a first look at some of the key fea- 
tures in the data. Plasma diagnostics of individual species will be 
presented in a later paper. 



2 DATA REDUCTION 

MCG-6-30-15 was observed by XMM-Newton on revolutions 
301, 302 and 303 (2001 July 31 - 2001 August 5) for 320 ks 
jFabian et alj|2002h . The observation data files (ODFs) were re- 
duced with version 5.3.3 of the XMM-Newton Science Analysis 
Software (SAS) using the standard processing chains and the cal- 
ibrations from 2002 December. The final few ks of data from each 
revolution showed strong background flaring, these periods were 
removed from the data extraction. The data were grouped so that 
at least 20 counts we re in each bin and were fitted to trial mod- 
els in XSPEC Vl 1.2 iArnaudll996T) . Data above 1.7 keV were ig- 
nored due to the known calibration uncertainties at the mirror Au 
edge. Unless otherwise stated errors quoted are for 90 per cent con- 
fidence level for one interesting parameter (A^ 2 = 2.7). The data 
are plotted in the observation frame throughout. 



3 SIMPLE ABSORPTION MODEL 

The fluxed RGS spectrum from all three revolutions is shown in 
Fig.Q Evident in the spectra are the drops in flux at 0.39, 0.54 and 
0.7 keV. The high spectral resolution of the RGS requires the in- 
clusion of individual resonance absorption lines when considering 
the effect of photoionized absorption. Specifically, at the red side 
of an absorption edge there is a series of resonant absorption lines 
that converge to the edge energy. These have the effect of reducing 
the observed edge energy when convolved with intrinsic Doppler 
broadening and the spectral resolution of the instrument. 

Models of the cross-section of OVII and OVIII near their K- 
edges were constructed using the con tinuum cross-section fitting 
formulae of IVerner & Yakovlevl 1 19951) . the resonance absorption 
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Figure 1. Combined and fluxed RGS data from all three satellite revolu- 
tions. 
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Figure 2. Detail of the predicted absorption near the OVII K-edge con- 
volved with the RGS resolution for an OVII column density of 10 18 cm -2 
and turbulent velocity of 100 kms -1 . (solid): Absorption including pho- 
toelectric edge and Is 2 — lsnp resonance lines, (dot-dash): Absorption 
including photoelectric edge only, (dash): Position of OVII K-edge. 

line data of lVemer et aljfl99 6) and the fi tting formulae for line ab- 
sorption series near absorption edges of Mewt] <1977l) . The trans- 
mission through a column density of 10 18 cm -2 of OVII, calcu- 
lated assuming a turbulent velocity (b parameter) of 100 kms -1 , 
is plotted in Fig. [2] This shows the effect of the resonance lines on 
the apparent position of the edge, which is clearly shifted to lower 
energies. Also present at 0.7 keV are the Fel L-shell absorption 
edges which have fin e structure in their cro ss-section at the L2,3 
edge as described bv lKortright & KimlfcOOOl) . 

Fig. [3] shows the revolution 301 data fitted to a power- 
law, Galactic absorption (TA^h = 4.06 x 10 20 cnT 2 ; 
Elvis. Wilk es. & LockmarJll989I) and absorption by OVII, OVIII 
and Fel across the energy range 0.6 - 1.7 keV. The best fit (\ 2 = 
9003.9 for 6625 degrees of freedom (dof)) requires a photon index 
of T = 2.36 ± 0.02 and column densities of log AW11 = 18.22 ± 
0.02, log/Vovin = 18.47 ± 0.02 and log 7V Pc i = 17.29 ± 0.02. 
Although the best fitting model gives a poor overall xt, it repro- 
duces the detailed shape of the drop at ~0.7 keV remarkably well 
as can be seen from Fig. [4] where the model is compared to the 
fluxed spectrum from all three orbits. The sharp drops at 0.7 and 
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Figure 3. Ratio of data for (filled circles) RGS 1 and (open circles) RGS2 
from orbit 301 to a model consisting of a power-law, Galactic absorption 
and OVII, OVIII and Fel absorption. 
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Figure 4. Details of fluxed spectrum for all three orbits compared to a 
model of OVII, OVIII and Fel absorption (solid) and Fel absorption only 
(dash-dot). 

0.713 keV are clearly fitted by the Fe I absorption profile. Also ev- 
ident are the resonant absorption lines of OVII. The Fel L3 edge is 
observed at 0.7023 ± 0.0007 keV which corresponds to an velocity 
of —320 ± 340 km s" 1 with respect to the rest fr ame of the source 
if the iron is in the metallic form described bv iKortright & Kiml 
(2000). The absolute error on the energy scale of measurements of 
the cross-section of iron oxides are required before it can be de- 
termined whether the observed iron is in the metallic or iron oxide 
form. The derived power-law photon index is higher than the pho- 
ton index of ~ 1.9 derived with the EPIC instrument above 2.5 keV 
(Fabian et al. 2002) suggesting the emission below 1-2 keV has a 
soft excess relative to an extrapolation of the best fitting continuum 
above 2 keV (see also lLee et all200ll) . The residuals around 1 keV 
are due to absorption lines of FeXIII-XXIV and NelX-X. 



4 PHOTOIONIZATION MODEL 

While the previous Section suggests that absorption by Fel, OVII 
and OVIII can explain the sharp drop in flux at 0.7 keV, it does not 
provide a self-consistent model of absorption from a photoionized 



absorber. In order to explore the properties of any warm absorp- 
tion in MCG-6-30- 15 and provide a model of this absorption the 
9 6beta5 version of the photoionization code CLOUDY was used 
to calculate the absorbed spectrum produ ced by the p assage of a 
continuum through a layer of gas iFerland et al. 1998). This ver- 
sion of CLOUDY was used since it includes th e atomic transitions of 
the ir on Unresolved Transition Array fUTA: lBehar. Sako. & Kahrl 
l200ll) which have been observed previously in Seyfert 1 X-ray 
spectra. Above 1 keV the incident continuum is taken to have the 
spectral shape inferred from fitting the EPIC and BeppoSAX data 
(Fabian et al. 2002). Below 1.0 keV the spectrum is assumed to be a 
power-law with V = 2.0 — 3.0 down to a break energy of 10 eV and 
T = —1.5 below this. The normalization of the incident spectrum 
is determined from the X-ray luminosity in the 2-10 keV range, 
L2-10, inferred from the EPIC and BeppoSAX data. Also included 
in the incident continuum is a blackbody component with a temper- 
ature of 2.5 xlO 5 K and bolometric luminosity of 6 x 10 43 ergs -1 . 
The gas is assumed to have a constant density and Solar elemental 
abundance with a turbulent velocity of 100 km s -1 . This turbulent 
velocity is used sinc e the lines appea r unresolved in the Chandra 
HETGS spectrum of lLee et alJfeOOll) . 

The code was used to generate a grid of incident and absorbed 
spectra as well as a list of absorption line optical depths at line 
centre where greater than 0.01. The parameters that were varied 
across the grid were column density of the gas (TVh = 10 20 — 
10 22 cm -2 ), inner radius of cloud from central source (R = 10 16 — 
10 19 cm), ionization parameter of the cloud (£ = 10 -5 — 10 3 ) and 
spectral index between 10 eV and 1 keV (r = 2.0 — 3.0). The 
density of the gas, tih, is set from R, £ and L2-10, through the 
formula: 

£ _ L2-W 

K ~ n H R 2 

The resulting transmitted continua (i.e. including only bound-free 
edges) were divided by the incident continua to produce a contin- 
uum absorption model. The absorption lines were then included 
using a Gaussian opacity profile with a b width of 100 km s~ . 

In using only the turbulent velocity as the line broadening 
mechanism the natural and thermal broadening of the lines are ne- 
glected. Thermal broadening is ignored since for the species of in- 
terest in a photoionized gas thermal broadening veloci ties are only 
10-20 per cent of the turbulent velocity iNicastro. Fiore. & Mattl 
ll99St) . lLeeeTail 1200 ll) show that various absorption lines from 
OVII and OVIII lie on the flat part of the curve of growth for col- 
umn densities of interest here, justifying our neglect of natural line 
broadening. 



5 SPECTRAL FITTING 

5.1 Dusty warm absorber model 

We now fit the warm absorber photoionization model to the data 
from all three revolutions. Galactic absorption, Fel absorption and 
multiple warm absorption zones are included to account for ab- 
sorption of the continuum. The absorption component is kept fixed 
between data from different revolutions while the parameters of the 
underlying continuum are allowed to vary. The continuum is taken 
to be a broken power-law with a break energy of 1 keV, spectral 
index of 2.2 above the break energy and a variable spectral index 
below the break energy. The use of a broken power-law is justi- 
fied since the variability suggests some break in the continuum at 
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Figure 5. Ratio of data from RGS 1 (filled circles) and RGS2 (open circles) 
to the best fitting DWA models. Only data for orbit 301 is shown here, (a) 
Four zone DWA model with underlying broken power-law. (b) Four zone 
DWA model with underlying broken power-law and excess broad emission 
at ~ 0.56 keV. 

~ 1 - 2 keV iFabian et alJl2002l : iFabian & VaugharJl2003l) . Fit- 
ting in the range 0.6-1.7 keV with a three zone warm absorber 
gives x 2 = 8412.6/6620 dof. The fit parameters for this model 
are given in Table Q (model 1). Three zones are required since 
two distinct veloc ity components are seen in the absorption lines 
ISako et aljEool and features from a broad range of ionization 
parameters are present. Extending the model over the whole RGS 
passband reveals that the predicted model flux is far higher than the 
observed data below 0.5 keV. The smooth nature of the discrepancy 
suggests additional cold absorption is needed, possibly in the form 
of absorption from the ISM within MCG-6-30-15. However, in- 
cluding additional cold absorption in our fit to the 0.32-1.7 keV en- 
ergy range produces a large OI edge that is not seen in the data. One 
possibility is that gas is partly ionized so most of the O is ionized to 
higher species reducing the depth of the OI edge. A four zone warm 
absorber produces a fit with x 2 = 16601.9/12053 dof, where one 
of the zones has a low ionization parameter log £ = —4.4 account- 
ing for the excess low energy absorption. Table Q (model 2) gives 
the fit parameters for this model. The residuals for this model are 
shown in Fig. [5^. As mentioned the lowest ionization parameter 
zone is required to cause the spectrum to turn down at low energies 
without too deep a neutral oxygen edge. The medium ionization 
parameter zone is required to fit the OVII edge and lines, the drop 
in flux at ~ 0.39 keV due to the CV edge and the Fe UTA at 0.7- 
0.8 keV which takes the form of a trough in the data. The higher 
parameter zones are required to fit a range of absorption line fea- 
tures around 1 keV as well as any OVIII features. The majority of 
the OVII column density is located in zone 2 (log£ ~ —0.5) while 
the majority of the OVIII column density is located in the highest 
ionization parameter zones 3 and 4 (log£ ~ 1.5). The best fitting 
DWA model contains less OVII and OVIII column density than the 
simple absorption model presented in section|3]since other absorp- 
tion components are present in the photoionization models, such as 
the Fe UTA, that can compensate for some of the OVII and OVIII 
column. Although residuals are present they take the form of broad 
features, whereas the region across the 0.7 keV drop is reasonably 
well fitted. 

We now add a broad Gaussian emission component at ~ 
0.6 keV to account for any excess emission in this region possi- 




0.66 0.66 0.7 0.73 0.74 

Observed Energy (keV) 



Figure 6. (a) Detail of residuals for the best fitting DWA model in the 0.65- 
0.75 keV region. The six data sets are shown simultaneously (one from each 
revolution and instrument), (b) Same as (a) except for the best fitting REL 
model. 

bly from reflection components. The best-fitting model produces a 
significant improvement in the fit with x 2 = 15977.3/12051 dof 
as can be seen in Fig.[5J) (see model 3 in TableQfor the fit param- 
eters of this model). The line has an equivalent width of ~ 120 eV 
and is very broad with a ~170 eV. The model fits the region 
around the 0.7 keV drop extremely well (see Fig. |S^). However, 
the broad Gaussian shape is not that of a single relativistic disc line 
and should be considered as only an approximation of the shape of 
the residuals. Adding the broad Gaussian emission line results in 
changes to the best fitting DWA component. The broad emission 
line has the effect of fitting the excess in the residuals in the range 
~0.5-0.9 keV. This produces an excess below ~0.4 keV and above 
~ 1 .35 keV. The DWA component compensates to fit these excesses 
by, at low energy, increasing the ionisation parameter of the lowest 
ionization parameter zone (zone 1) which results in absorption that 
rolls over less at low energies and is flatter, and at high energies 
by decreasing the column density of the high ionization parameter 
zones which results in less absorption above ~ 1.35 keV flattening 
the absorption. 

5.2 Relativistic emission line model 

iBranduardi-Ravmont etailfeOOll) and lSako etal]j2003l) argue that 
the drop in flux at 0.7 keV is due to the blue wing of a rela- 
tivistically blurred OVIII Lya emission line. In this interpretation 
the warm absorber is responsible for only narrow absorption fea- 
tures rather than any broad ionization edges. In this section we fit 
this relativistic emission line (REL) m odel to data from all three 
revolu tions. In the original model of IBranduardi-Ravmont et ail 
(2001) the model consists of a power-law and emission lines due 
to CVI Lya (0.368 keV), NVII Lya (0.500 keV), and OVIII Lya 
(0.654 keV) emitted by the inner-part of an accretion disc around 

a rotating black hole. The relativistic blurring is achieved with the 

I 1 1 I 

LAOR model <Laorlll9 91), Free parameters in the model are in- 
ner and outer radius of the disc (7?i n and 7? ut), disc inclination 
(i), emissivity e as a function of radius R which is parametrized 
as a power-law with index q (i.e. e cx R~ q ). Galactic absorption 
is also included. The best fitting model produces a fit with x 2 — 
17457.8/12051 dof (see Fig.|7Ji). The fit parameters of this model 
are given in Table|2|(model 1). While the level of the continuum is 
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Table 1. Results of fits of the DWA model to all three revolutions of RGS data. Line energy of the broad Gaussian emission line (-Eg) is given in units of keV. 
Width (cr g ) and equivalent width (EW S ) of the broad Gaussian emission line are given in eV. Column density (JVx) is given in units of log cm -2 . Velocity 
(v) is given in kms -1 . 



model 


zone 


1 


2 


3 


l"i 




2.72 ± 0.03 


2.91 ± 0.02 


2.68 ± 0.04 


Ea 








0.54 ± 0.03 


CT g 








176 ± 15 


EW g 




- 




121 ± 12 






17.28 ± 0.03 


17.22 ±0.02 


17.25 ±0.03 


Ytt 
1 v 11 


1 


21 37 + 04 

_±.'Ji _L_ W»V7^t 


20 86 ± 01 


21.12 ± 0.04 


log 5 


1 


-0.47 ± 0.03 


-4.42 ±0.23 


-1.92 ±0.02 


V 


1 


0/ 


0^ 


0^ 




2 


21.72 ± 0.05 


21.44 ± 0.03 


21.44 ± 0.03 


log? 


2 


1.43 ±0.03 


-0.52 ±0.02 


-0.36 ±0.03 


V 


2 


0/ 


0/ 


0/ 


Nu 


3 


22.00-f 


22.00-f 


21.73 ±0.10 


log? 


3 


2.00 ± 0.04 


1.69 ±0.03 


1.73 ± 0.06 


V 


3 


-2200 ± 100 






N H 


4 




22.00-f 


21.29 ±0.08 


log? 


4 




1.78 ± 0.03 


1.39 ±0.05 


V 


4 




-1730 ± 50 


-1910 ± 140 


Novn 




17.81 


18.11 


17.95 


NoviII 




18.14 


18.20 


17.97 


X 2 /dof 




1.271 


1.377 


1.243 






(8412.6/6620) 


(16601.9/12053) 


(14977.3/ 12051) 



' -fixed. 



well reproduced, the region across the 0.7 keV drop is poorly fitted 
since no Fel L-shell absorption is included. Adding Fel to account 
for the L2.3 edge clearly seen in the data (AW = 10 17,2 cm~ 2 ) 
produces a fit with x 2 = 17534.6/12051 dof. See Table|2l(model 
2) for fit parameters. The effects of warm absorption are now in- 
cluded by adding absorption by OVI-OVIII, FeXVII-XXIV and 
NelX-X in two velocity components, one with an outflow veloc- 
ity of ~ 2000 kms -1 an d the other at rest in the frame of the 
source as was observed by Sako et al. 1 2003). Also included is a 
broad Gaussian absorption trough at ~0.75 keV to represent any 
Fe UTA. Allowing the columns of these species and Fel to be free 
parameters improves the fit (\ 2 = 17579.4/12022; see Fig.[7J)). 
Table[2|(model 3) gives the fit parameters for this model. The OVIII 
column density of 10 18,4 cm -2 produces an edge at 0.871 keV with 
t ~ 0.24. The best-fitting column of Fel required is insufficient to 
account for the fine structure seen at 0.7 keV (see Fig.|6j)) and the 
spectral index of the power-law is ~1.5, muc h lower than the ~2.0 
index found for the EPIC data I Fabian et al. 2002). The equivalent 
widths of the three emission lines are ~34 eV for CVI Lya, ~38 eV 
for NVII Lya and ~220 eV for OVIII Lya. 

Under the REL model of ISako et alJ <2003l) the 1.0-1.7 keV 
region consists of a power-law modified by narrow absorption fea- 
tures. Here we fit this region to a power-law producing a best-fitting 
spectral index of F = 0.93 ± 0.04. Even with a column of Fel of 
10 17,2 cm~ 2 the best-fitting spectral index is 1.10 ± 0.04. In order 
to match the continuum here with the softer power-law found at 
higher energies the continuum is required to break at around 2 keV 
from a hard power-law below the break to a softer one above. No 
mechanism is known that can create such a break. 
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Figure 7. Ratio of data from RGS 1 (filled circles) and RGS2 (open circles) 
to the best fitting REL models. Only data for orbit 301 is shown here, (a) 
Relativistically blurred emission lines and power-law (b) Relativistically 
blurred emission lines and power-law with warm absorption. 



5.3 Combined DWA/REL model 

Finally the DWA model was combined with the REL model and 
fitted to the data. It was found that the best fitting model obtained 
from the fit was strongly dependant on the starting position of the 
fit. If the REL model with zero line normalisation was added to the 
best fitting DWA model (model 2) then only very small emission 
lines were required. When the DWA model was added to the best 
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Table 2. Results of fits of the REL model to all three revolutions of RGS data. Equivalent widths of the emission lines (EW) are given in eV, disc radii (ijj 
and -Rout) are given in units of r g , inclination (i) is given in degrees and column density {Nh) is given in units of log cm -2 . 



model 


1 


2 


3 


r 


1.46 ±0.02 


1.75 ±0.03 


1.50 ±0.04 


EWcvi 


35.2 ± 12.7 


27.6 ± 5.2 


34.1 ± 10.0 




KG K -|_ q K 


i k 2 + 3 8 

lU.i _1_ ■ ) . O 


QQ A J- 1 K K 


-C/ vv OVIII 


417.3 ± 7 


167.7 ± 4.2 


21 s ^ -J- q ■} 


q 


3.81 ± 0.02 


4.68 ± 0.07 


3.79 ± 0.08 


R 


9 49 + n (11 


i on + n 27 


2 11 + n 02 


Rout 


185 ± 85 


13.69 ±3.60 


11.2 ±0.3 


i 


37.3 ±0.15 


45.1 ±0.18 


40.0 ±0.2 


iV H (FeI) 




I12 f 


16.78 ±0.03 


Af H (OVII) 






17.32 ±0.08 


N H (OVIII) 






18.39 ±0.03 


X 2 /dof 


1.449 


1.455 


1.310 




(17457.8/12051) 


(17534.6/12051) 


(15749.4/12022) 



f -fixed. 



fitting REL model (model 1) then very little warm absorption was 
required. It would appear that two separate solutions exist; one with 
a large DWA and small REL and the other with small DWA and 
large REL. 

5.4 DWA/REL fitting summary 

The REL model by itself is a poor fit and requires the addition of 
a DWA. Even then some parameters are unphysical. The simple 
DWA model, acting on a power-law continuum, probably requires 
structure in the soft X-ray continuum. The likely best-fitting model 
is a DWA acting on a continuum with some emission structures. 



6 EMPIRICAL ABSORPTION MODEL 

The principal difficulty in fitting the soft X-ray spectra of absorbed 
Seyfert 1 AGN is that the shape and level of the continuum is not 
known, partly because broad emission lines can mimic absorption 
edges. In this section we attempt to circumvent this problem using 
the known spectral variability of MCG-6-30-15. 

The flux in the Fe Ka line at ~6.4 keV is observed 
to va ry little w hile the flux in the continuum is highly vari- 
able <Lee et alJl2000t Ivaughan & EdelsoiJl200ll; Lee et alj|2002t 
IShih. Iwasawa, & Fabian 2002; Fabian et al. 2002). If a spectrum 
extracted from a period of low continuum flux is subtracted from 
a spectrum extracted from a period of higher flux, the resulting 
spectrum (the difference spe ctrum) can be adequately fit by just 
a power-law above ~2 keV iFabian et alj l2002). It would there- 
fore appear that above ~2 keV the emission consists of a rapidly 
varia ble power-law with constant Fe Ka emission. On this basis 
iFabian & Vaughanl (2003 ) introduce a model for the emission con- 
sisting of a rapidly variable Power-Law Component (PLC) and an 
almost constant Reflection Dominated Component (RDC). Since 
the RDC changes little, subtracting a low flux spectrum from a high 
flux spectrum effectively removes it. 

Data from revolution 303 were used to examine this variabil- 
ity for the RGS data since the flux level varied considerably during 
this period (see Vaughan et al. 2003 for light curves). Two pn spec- 
tra were extracted, one from the first 80 ks when the mean flux 



level was high and one from the last 40 ks when the mean flux 
level was low. The low flux spectrum was subtracted from the high 
flux spectrum and the resulting difference spectrum was found to 
be adequately fitted by a power-law with T ~ 2.2 in the range 
3-10 keV. The ratio of data to best fitting power-law is shown in 
Fig. UK Assuming that the variability model applies to any reflec- 
tion components below 3 keV, the underlying continuum in the pn 
difference spectrum is just a power-law and any deviations that are 
seen in Fig.|8ji are due to absorption. These deviations from the 
high energy power-law therefore represent the absorption function 
of the source, including both intrinsic and Galactic absorption. 

Next RGS spectra for the first 80 ks and last 40 ks of revolu- 
tion 303 were extracted. The RGS spectra from the low flux period 
was subtracted from the RGS spectra from the high flux period for 
each instrument to leave the PLC modified by absorption. These 
RGS difference spectra were then divided by a power-law, with the 
spectral index and normalization fixed to the values derived from 
the pn difference spectrum, to produce the RGS absorption func- 
tions in Fig.|8j) and c. The OI edge at ~0.54 keV and the Fel/OVII 
edge at ~0.7 keV can clearly be seen. This absorption function is 
compared to the best-fitting DWA absorption model from section 
15. ll in Fig.|9jt where the best-fitting DWA model is seen to be more 
absorbed than the absorption function inferred from the pn differ- 
ence spectrum. This would result if the varying PLC is a broken 
power-law with a break at ~ 1 — 2 keV and a softer spectrum below 
this break than that inferred from the pn difference spectrum in the 
3-10 keV range. Fig. |9j> compares the empirical absorption func- 
tion obtained by assuming the PLC to be a broken power-law with 
slopes Ti = 2.66, T 2 = 2.29 and break energy E hl = 1.75 keV, 
with the best-fitting DWA model. Good agreement is found. Also 
if the best fitting DWA model is applied to the pn difference spec- 
trum a broken power-law for the underlying continuum produces a 
good fit to the data. Fig. llOl shows the different model components 
present in the spectra. 

We proceed by assuming that the absorption function for the 
soft X-ray spectrum is the pattern of deviations in the RGS dif- 
ference spectrum from a broken power-law. This will be the case 
if the variable part of the spectrum is approximately a broken 
power-law continuum. Specifically this means that any additional 
emission features do not vary with the continuum, which would 
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Figure 8. (a) Data/model ratio plot from fitting the pn difference spectrum 
to a power-law in the range 3-10 keV, excluding Galactic absorption, (b) 
Data/model ratio plot for the RGS 1 difference spectrum compared to the 
model from (a), (c) Same as (b) except for the RGS2 difference spectrum. 
The gaps in the RGS data are due to nonfunctional chips in the detectors. 
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Figure 9. Empirical absorption function derived from RGS 1 (filled circles) 
and RGS2 (open circles) compared to the best fitting DWA model absorp- 
tion model (solid line), (a) Absorption function derived from the pn differ- 
ence spectrum, (b) Absorption function derived from a broken power-law. 

mean they contribute to the reflection dominated component of 
iFabian & Vaughanl J2003 ) . The assumption is supported by the fact 
that the RMS variability spectr um shows no shar p emission or ab- 
sorption features below 1 keV (Fabian et al. 2002) and is well rep- 
resented by the above model. The spectral indices and break energy 
of the broken power-law are assumed to be fixed while the normal- 
ization varies. Such a continuum is observed in many other Seyfert 
1 X-ray spectra. The RGS difference spectra are divided by the bro- 
ken power-law to form the RGS absorption function. Now that the 
absorption function has been derived for the soft X-ray band, its ef- 
fects can be removed from the RGS data for the whole observation 
and the RDC in the soft-band can be examined. RGS data for all 
three revolutions were divided by the absorption function to cor- 
rect for both Galactic and warm absorption. The spectra for each 
revolution were then fitted to a power-law and the ratio of data to 
best fitting model for each revolution combined to produce Fig. II II 
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Figure 10. Model components present in spectra, (a) Power-Law Compo- 
nent (PLC), (b) Reflection Dominated Component (RDC) showing strong 
Fe Ka emission and weaker soft X-ray emission lines, (c) Dusty Warm Ab- 
sorber (DWA) showing strong OVII and OVIII edges, (d) Combination of 
all model components. N.B. This diagram is purely schematic and does not 
represent the actual model components. 

The photon index of the best fitting power-law for the absorption 
corrected spectra was ~2.8, which is steeper than that derived for 
the high energy spectral index of the PLC of ~2.3, re-enforcing the 
idea that there may be a soft excess below 2 keV, similar to t he vari - 
able two component power-law in Ton SI 80 <VaughanetaiT2 002). 
The residuals also show an excess of emission below 0.4 keV and 
a trough at ~1 keV. Any features around the OVII (0.5740 keV) 
or OVIII (0.6536 keV) emission lines are constrained to be less 
than ~20 per cent of the continuum. Provided the PLC continuum 
is a power-law in the RGS passband, the shape of the residuals in 
Fig. II H are independent of the spectral index and normalization as- 
sumed for that power-law. Note the resemblance to the residuals 
found ear tier (Fig.EJi). 

The whole emission complex between 0.45 and 0.9 keV was 
earlier modelled as a broad Gaussian (section |5~TV We now fit it 
in the absorption corrected spectra with a power-law and emission 
lines centred on the OVII Is 2 — ls2p and OVIII Lya line energies. 
Fits with both intermediate width Gaussian and relativistically- 
blurred lines produce upper limits on the equivalent width of OVII 
and OVIII emission of ~ 30 eV. The whole complex could be due 
to OVII, OVIII, FeL and other em ission expected in the reflection 
spectrum iBallantvne et alJl2002l) together with further absorption 
structure. Detailed modelling of this emission will be carried out in 
later work. 



7 DISCUSSION 

The high signal-to-noise of our 320 ks observation has allowed the 
unambiguous detection of Fel L-shell absorption through the fine 
struct ure at the L2.3 edg e. This has been observed before in th e 
ISM IPaerels et alfeoOlh and in Cygnus X-l ISchulz et alJl2002l) . 
It combines with the resonance and edge absorption from OVII and 
OVIII to explain the drop in flux at 0.7 keV (Lee et af]l200ll) and 
remove the ~ 16,000 kms~ discrepancy in position between the 
observed drop and that expected from the OVII K-shell edge re- 
ported in Branduardi-Rav mont et alJ (2001 ). 

In comparing the DWA and REL models to the data, both pro- 
duce similar values of x 2 > although all models are formally un- 
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Figure 11. Combined ratio of the absorption corrected data from RGS1 
(filled circles) and RGS2 (open circles) to the best-fitting power-law for all 
three orbits. 



tinuum has some complex emission structures in the 0.4-0.9 keV 
band of amplitude consistent with theoretical predictions for an ion- 
ized reflector. 
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acceptable. However the normalizations of the three broad emis- 
sion lines in the REL model are free parameters, so the level of 
the continuum in each section of the spectra can be easily fit. 
The equivalent widths of these lines are also much larger than has 
been obtained with self-consistent models of ionized disc reflection 
I Balla ntvne etai]|2002h . As has been noted in Section|5| the best- 
fitting REL model requires insufficient Fel to fit properly the fine 
structure at 0.7 keV and produces an unrealistically flat spectral in- 
dex for the underlying power-law. The DWA model in contrast uses 
only a broken power-law for the underlying continuum. Addition 
of a reflection component, with emission features around the 10- 
20 per cent level fFie. ll II . may well exp lain some of the rem aining 
residuals (see for example the model in Fabian & Vaughan 2003h . 
The DWA model also assumes Solar abundances and a turbulent 
velocity of 100 kms -1 . Relaxing these two assumptions should 
improve the fit further. These improvements will be explored in a 
later paper. 

In Section [5] the pn and RGS difference spectra are used to 
determine the presence of additional emission components. This 
method reveals a steep increase in flux below 0.4 keV and a trough 
at 1 keV. This suggests that an extra emission component is present 
below 0.4 keV which, combined with increased low ionization pa- 
rameter warm absorption, could explain the low energy residuals 
from the DWA model. This extra emission component may be the 
high energy tail of the quasi-blackbody emission emitted from the 
accretion disc. These methods assume the absorption present is 
non-varying, whereas an anti-correlation of the optical depth of an 
absorption feature with ionizing flux would produce a drop in the 
absorption corrected spe ctra where that abs orption was. Such vari- 
ability was observed bv lOtani et alj (1996) for absorption around 
1 keV. The residuals in the DWA fit at 1 keV may be due to a warm 
absorber zone with this sort of variability which would explain the 
trough seen at ~ 1 keV in the absorption corrected spectra. The 
fact that no large peak in the absorption spectra is observed in the 
0.6-0.7 keV region suggests that any OVIII emission, which does 
not vary exactly with the continuum, is at a level not greater than 
~ 10-20 per cent of the continuum. This is much smaller than in 
the basic REL model. 

In summary, we have found that detailed absorption modelling 
of the average RGS spectra agrees with an analysis based on vari- 
ability. Both indicate the dominant presence of dusty warm ab- 
sorber. Residuals from the modelling show that the underlying con- 
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